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Signal Processing Techniques for Clutter Filtering and Windshear 

Detection 

E. G. Baxa, Jr., Clemson University 
ABSTRACT 

It has been argued that the windshear hazard factor is a 
sufficient statistic for detecting hazardous windshear conditions. 
The hazard factor is computed by estimating the spatial gradient 
of windspeed across the radar sector of coverage. With the 
airborne Doppler radar, one approach is to use estimates of 
windspeed within each range resolution cell as a basis for 
estimating this spatial gradient. Currently, research is directed 
at understanding how to obtain the best possible estimate of 
windspeed conditions within a range cell. Conventional pulse-pair 
processing obtains mean estimates of windspeed. The presence 
strong ground clutter in a low altitude airborne radar return can 
significantly bias these mean estimates. One thrust of this 
effort has involved use of adaptive clutter rejection filters 
based upon auto - regressive modelling of the ground clutter 
returns. This offers the potential for using very simple finite 
impulse response digital filters to eliminate highly specular 
ground clutter returns. For situations where the weather return 
is quite low, e.g., the "dry" microburst, clutter rejection 
filtering can reduce the weather return signal levels to the 
extent that the variance of the mean estimates is quite large. 
Research is involved with using mode estimates, i.e., estimates of 
the most probable windspeed, in each range cell in determining the 
hazard factor. An extended Prony algorithm is discussed. It is 
based upon modelling the radar return as a time series and appears 
to offer potential for improving hazard factor estimates in the 
presence of strong clutter returns. 
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FINAL APPROACH SITUATION 


dry microburst 



Simulated final approach situation with 
A/C on 3 degree glideslope and radar 
antenna elevated 2 degrees . Dry 
microburst in front of Denver runway 26R. 

Ground clutter return is based upon SAR 
data taken at Denver Stapleton airport. 

Signal to clutter ratios are on the order 
of 0 dB in the range cells in which the 
microburst is present. 
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DRY MICROBURST 


10th Order AR Spectrum 
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Auto - regressive model determined spectrum 
in each of the fifty range cells with 
the simulated "dry 11 microburst without 
any clutter present. Signal* to -noise 
ratios in the range cells with the 
microburst varies from 0 to 30 dB. 

Note: zero windspeed corresponds to zero 
Doppler relative to the ground speed 
of the aircraft. Positive windspeed 
corresponds to winds toward the aircraft 
and negative is away from the aircraft. 
Range cells are 150 m. 


DOPPLER MEAN ESTIMATES 
weather only - mean estimates 

10th order 
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Mean windspeed estimates considering 
simulated "dry" microburst without 
ground clutter. Five different mean 
estimates are used: 

1. pulse-pair computed in the time domain 

2. pulse-pair computed in the frequency 
domain using an AR spectrum estimate 

3. Fourier domain mean estimate 

4 . AR spectrum domain mean estimate 

5. First order AR model pole estimate 

Note: The microburst appears in range 
cells 20-33 (approximately). Some 
estimates of mean have been edited 
to zero outside this range based upon 
estimated signal to noise ratio in return 
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ar 

Width estimates for the situation in the 
previous slide. Four different width 
estimators have been used: 

1. pulse-pair width computed in the time 
domain . 

2. pulse -pair width computed in the AR 
spectrum frequency domain. 

3. AR spectrum standard deviation 

4. First order AR model coefficient 

Note: Some width estimates for range cells 
outside those containing the microburst 
have been edited to zero because of low 
signal to noise ratio estimates. 
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DOPPLER MODE ESTIMATES 
weather only 

2nd order extended prony mode estimates 

2 largest modes 
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Spectrum mode estimates using an extended 
PRONY algorithm based upon a second 
order AR model of the data. Outliers 
are caused by insufficient model order. 
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DRY MICROBURST IN CLUTTER 


10th Order AR Spectrum 
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AR model determined spectrum in each of 
the fifty range cells with the "dry" 
microburst and ground clutter present 
in the return. No clutter rejection 
filtering is used. Ground clutter in 
the range cells 40-50 in the negative 
Doppler region is associated with 
an interstate highway included in 
the simulation. Signal to clutter 
ratios are on the order of 0 dB. 

Note that the microburst can still be 
identified. 




Notes 


Mean estimates without clutter rejection 
filtering for the situation depicted 
in the previous slide. The same five 
estimators used proviously are included. 
Again some of the mean estimates have 
been edited to zero based upon signal 
to noise ratio estimates of the return. 

Note that the clutter biases the mean 
estimates in the range cells 20-33 
so that the presence of the mocroburst 
is no longer evident. 
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DOPPLER MODE ESTIMATES WITH CLUTTER 
weather + clutter 

10th order extended prony mode estimates 

2 largest modes 
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Spectrum mode estimates using an extended 
PRONY algorithm based upon a tenth 
order AR model of the data. Only the 
two strongest modes within each range 
cell are retained. Outliers are 
caused by the presence of discrete 
clutter (e.g. interstate highway) 

Note that the microburst spectrum modes 
are clearly identifiable even though 
no clutter rejection filtering has 
been done. 





POST-FILTERED MEAN ESTIMATES 
optimally filtered clutter+ weather - mean estimates 

1<Xh order 
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Mean estimates with optimum clutter rejectioi 
filtering. The same five mean estimators 
used previously are compared. Again 
some of the mean estimates have been 
edited to zero because of low signal 
to noise ratio estimates. 


Note that the microburst can be clearly 
identified. 





PHASE NOISE EFFECTS ON CLUTTER FILTERING 
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Notes 


Radar system pulse- to -pulse phase jitter 
is analyzed in the presence of the 
low signal to clutter ratio situation. 

Here an ideal notch filter centered at 
zero Doppler with a stopband width 
of 150 Hz and 70 dB stopband attenuation 
is analyzed. The prefiltered signal 
to clutter ratio is held to -30 dB 
and the weather mean Doppler is varied 
from 100 to 250 Hz. As the phase jitter 
noise is increased the clutter spectrum 
is spread to the point that the rejection 
filter will not provide enough signal 
to clutter ratio gain for reliable 
pulse pair processing. 
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SUMMARY 


Characterization of winds peed within a radar rang# resolution 
cell can be severely limited by ground clutter returns 


Low level weather returns will present the greyest challenge 
In Hazard detection 


Signal processing needs include a variety of algorithms 
and may require super-computer processing toads for 
real time implementation 
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Characterization of ground clutter returns 
in initial flight tests will be of 
paramount importance. 

A suite of signal processing algorithms 
will be needed to improve confidence 
in hazard detection 


Airborne radar will be important for 
hazard detection but should be 
integrated with other sensor types, 
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Estimation of Radial WindSpeed: 


* From the I and Q data, the mean radial windspeed 
is determined using Covariance and Spectral 
domain approaches. 

^ Here we study the pcrlormance of each of these 
techniques under varying signal to noise ratio. 


Covarian ce Method: 

If R(t) is the covariance function of the received 
sequence then the mean Doppler frequency f d can 
be estimated by 


fdS 1 


2nTr 


A , / lm -(R(Tr\ 

A n< Re(R(Tr) 


The mean radial wind speed is then obtained as 
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Spectral Estimation Methods 

If S(f) is the spectral density ot the sequence 
then fd can be estimated by using 

N/2 

I fiS(fi)W(fi) 

fd=i^lZ2 

N/2 

I S(fi)WCfi) 

i=-N/2 

where w <0 is the weighting function introduced to 
suppress the stationary ground clutter which is 
centered around zero Doppler frequency. 

The spectral density S(f) is determined using 
following methods. 


(1) Periodogram Method 

( 2 ) Forward-Backward Linear Prediction 
Method 

(3) Eigenvector method 


(4) MUSIC Method 


The performance of these methods when the 
signal to noise ratio is varied between 10 dB to -5 
dB.is studied 

From these results it may be concluded that 
Covariance method under severe SNR performs 
better than other methods. 
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PULSE PAIR PROCESSOR (SNR = -5DQ) 
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EIGENVECTOR METHOD (SNR — — 5 DR) 
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MUSIC METHOD (SNR=1QDU) 
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PERIODOGHA M METHOD ( SNR — INFINITY) 
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PERIODOGRAM M ET HOD (SNR=-5DB) 
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